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A  class  of  data  windows  is  proposed  for  use  in  the*  estimation  of  the  power 
spectra  cf  stationary  stochas»tic  processes.  These  windows,  which  fire  p.onoral i .'ra¬ 
tions  of  the  standard  Far/.c»»  filter,  me  constructed  by  the  use  of  appropriately 
normalized  B-splinos.  It  is  demonstrated  how  this  window  c  1  a.;s  nay  be  ce.puto- 
tiorally  impl »*’.»ent ed  using  c  last  Fourirr  Jrcnsi'orn  alyorithm,  Iho  efficiency  cf 
the  resulting  ptocedu.ro  is  genera!  lv  a  significant  improvement  over  the  stale  of 
the  art  with  little  additional  computer  tine*  required. 


i.  iNTRomrcrroN 


TM 5  psprr  f «?  concerned  with  the  prohloi,  of  the  result  of  unin^  un.se  cfU  windows  is  a  rig- 

digitally  estimciing  the  pever  spectrum  of  a  wide-  nlficant  fnpro/e-ient  ir  estimates  of  power  spectis 
sense  stationary,  ereodic,  Stochastic  process  without  a  noticeable  increase  in  the  commitat » c»nal 

from  a  sample  function  of  finite  length.  Not  effort.  We  have  chosen  to  u>c  continuous  notation 

long  after  Cooley  and  Tuhey  [5]  introduced  the  throughout;  conversion  to  discrete  formulae  where 

fast  Fourier  transform  algorithm  (FFT>,  its  appli-  appropriate  is  straight  forward, 
cation  to  the  estimation  of  power  spectra  was  dis-  2  PREJ IM  IN  ARIFS 

cussed  by  Binghara,  Godfrey,  and  Tukey  f 2 1  and  by 

Welch  [18],  They  suggested  a  method  based  on  com-  We  first  outline  some  basic  notation  and  re- 

plex  demodulation  and  pointed  out  the  computat ional  suits  which  are  assumed.  Here  we  follow  Stein  and 


effort.  We  have  chosen  to  me  continuous  notation 
throughout;  conversion  to  discrete  formulae  where 
appropriate  is  straight  forward. 

2.  PRELIMINARIES 

We  first  outline  some  basic  notation  and  re¬ 


speed  advantages  which  can  be  obtained  by  using 
the  FFT  for  computing  Fourier  periodograms .  Welch 
proposed  a  direct  method  which  has  become  the 
standard  numerical  approach  to  the  problem  of 
powe**  spectral  estimation.  This  procedure  has 
three  parts:  subdividing  the  available  data,  using 
^  ^  a  data  window  t<’  compute  a  smoothed  spectral  csti- 
w  mate  ^or  segment,  and  averaging  these  smoothed 

spectral  estimates.  The  present  paper  proposes 
{  ^  a  class  of  da fa  windows  derived  from  the  B -spline 
basis  functions  of  Schoenberg  (141,  Krppi rlca 1 1 v, 

,  . 


Weiss  (16).  I.1  denotes  the  Banach  space  of 

Lcbesgue  measurable  functions  which  arc  (Lcbesgv.c) 
integrable  over  the  real  line  H  .  For  f^L1  ,  tlie 
Fourier  transform  of  f  is  the  function  defined 
by 

ho  =  r„  f(t)C'2,Tiadt 

for  all  Z £ P  where  l  -  /-I  .  The  notation 
(...)*  denotes  the  Fourier  transform  of  (...). 

For  f,g€U  ,  the  convolution  h  -  f*p  Is  defined 
as  the  function  h  (  l.1  given  by 
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b(’>  «  f(*-t)g(t)dt 

for  .ill  ?6P.  Tf  f,grL  ,  then  «  fg. 

Similarly,  (fg)*  -  f*g  provided  the  functions  in¬ 
volved  are  such  that  all  operations  make  sense. 

For  a  >0,  fi4  represents  a  dilation  operator, 
i.o.,  \f(t)  *  f(at).  Whenever  f ,  gh  Ll,  it 

is  clear  that 

C.f*6.s)(n  =  l  \((f*g)(?)l 

and  (S,f)*(?)  =  i«.v  f(t)  . 

We  shall  consider  a  sample  function  x(t) 

T  T 

on  a  finite  interval,  '  2  -  C  -  2  unless 
otherwise  stated,  of  a  stochastic  process 
(X(t) : -« <  t < m) .  We  assume  this  underlying 
stochastic  process  to  be  wide-3ense  stationary 
(f)oob(8J)  and  ergodic.  For  simplicity  we  take 
the  mean  value  to  be  zero.  Let  X(c)  have  the 
spectral  density  function  P(f).  Hie  terminology 
related  to  spectral  analysis  used  in  this  paper  is 
generally  that  of  Blackman  and  Titkey  *[3). 

It  is  apparent  that  high  resolution  and 
stability  (snail  bias  and  variance)  are  desired 
qualities  of  power  spectral  estimates.  However, 
a  compromise  ween  the  two  must  he  made.  fn 
order  to  obtain  a  stable  spectral  estimate  using 
the  classical  indirect  method,  an  appropriate 
function  of  lag  must  modi fy  the  sample  autcco- 
variance  function.  Fundanenta 1 ly  the  indirect 
method  estimates  Pff)  via 

i’.  (f)  -  f  C(t)t.(  r)e'2rlfT  c!t 

t  .  M  .  w  _  T 

whore  l.(T)  is  the  finitely  supported  lag  window 
and  is  the  sample  aulocovar i ance  function 

t  ,(T-|t|)/2  _ 

C(T)  r -cr-lrh/2  x(t  +  )at- 

Usual  ly  the  support  of  I.(r)  Ls  much  less  than 
the  entire  interval  -T<  r<  T  . 

The  arguments  for  using  lag  windows  (e.g,, 
see  Jenkins  and  Watts  f 10))  ran  he  adapted  to 
show  that  a  p:oct*dure  which  utilizes  the  KFT  to 
estimate  power  spectra  directly  should  incorporate 
a  data  window.  That  is,  for  n  direct  procedure 
to  have  a  proper  balance  between  resolution  and 
stability,  the  data  itself  must  be  modified  by  a 
mltablr  function  of  time.  Thus ,  basically  the 


at  discrete  frequencies  where  W(t)  is  the 
finitcLv  supported  data  window. 

The  end  result  of  either  of  these  modifica¬ 
tions  is  that  we  in  fact  estimate  smoothed  values 
of  P(f).  More  specifically,  if  we  let  P„(f) 
represent  the  estimate  computed  by  either  the  in¬ 
direct  or  the  direct  approach,  we  have 

avg[  Pe  (f  )j 55  Q(f)*P(f) 

where  the  average  may  be  taken  either  over  the 
ensemble  or  along  time  [3J.  <)(f)  is  called  the 

spectral  window  corresponding  to  the  window  used 
in  the  computation. 

The  direct  and  indirect  approaches  to  esti¬ 
mating  the  power  spectrum  of  a  stochastic  process 
are  not  fur.da,  venta  1  ly  different.  Blackman  and 
Tukev  [31  have  derived  the  ccv.oct  ic"»  between  a 
data  window  Wft)  and  its  j  i  /a  lent "  lap,  win¬ 
dow  L(~).  'his  relationship  turns  out  to  he  the 
correlation  integral 

l.(-t  -  Wfnw(t-‘rVit  . 

By  "equivalent"  we  mean  that  rh*_.  appropriate  use 
of  the  respective  'Vquival  nt”  windows  yields 
power  spectral  estimates  which  jre  equal  in  .x- 
poctation.  Note  that  if  the  data  window  is  an 
even  function  of  time,  then  the  equivalent  lag 
window  is  proportional  to  the  convolution  of  the 
data  window  with  itself. 

At  this  point  it  shou  Id  he  clear  that  in  the 
indirect  method,  the  spectral  window  i}(0  cor¬ 
responding  to  the  lag  window  I. (-I  is  simple  the 
Fourier  transform  of  I-frl.  For  the  direci  method 
it  tii c * m  follows  that  the  spectra!  window  Q(f) 
corresponding  Lo  the  data  window  W(t  )  is  give*' 
bv  Q  (  f )  =  where  !  (  f  )  is  defined  to 

be  the  Fourier  transform  ot  'v(l).  J(f>  is 
called  the  frequency  window  eo» respond i  iv.  to 
W(t).  We  again  refer  the  reader  to  Blackman  and 

Tukey  (3). 
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3.  srur.r  data  ;:»N*novK 
i.l  motivation 

lr  light  of  the-  rather  simple  i olat ionshi p 
between  a  data  window  and  its  equivalent  Jar  win¬ 
dow,  one  should  not  ignore  the  i*>:tens ive  eifcit 
which  has  boon  (ml  into  tlu»  construction  of  lag 
windows  hv  numerous  ant  hots  nor  the  proper tits  of 
Lhose  windows.  See  for  example  the  cor;?i  Jr.  ?  :  r  n 
given  hy  Jenkins  f  9)  or  by  larzen  [12,13). 

For  instance,  recall  the  continuity  clays  to 
which  these  functions  belong.  In  particular  we* 
observe  that  one  of  the  lag  windows  which  has  been 
proposed  by  Parzen  f 13]  is  proportional  to  a  di¬ 
lated,  fourth  order  B-splino  basis  function.  One 
of  the  two  data  windows  suggested  by  Welch  118] 
has  the. shape  of  a  second  ord«r  B^-spline  basis 
function;  he  recognized  that  his  proposed  window 
yields  the  Parzen  spectral  window.  Henceforth 
we  will  refer  to  this  window  which  is  given  by 


as  Inc  Parzen  data  window,  even  though  it  is  also 
a  member  of  the  family  of  spline  windows  proposed 
in  this  paper.  Today,  Welch's  procedure  (de¬ 
scribed  In  the  introduction)  i ncorporat i ng  the 
Parzen  window  is  undoubtedly  the  most  widely  used 
method  of  direct  power  spectral  estimation. 
Bingham  et  al  (2)  proposed  the  following  data 
window  which  is  continuous  through  the  first  dc- 
rivat i ve: 


It  I <  .4T 


W(t)  =  < 


.At  <  It  I  <  .st 


^  0  ,  otherwise  . 

Wc  have  normalized  this  window  properly  although 
this  was  not  done  in  ( 2 ] . 

From  another  point  of  view,  we  note  that 
under  the  Fourier  transform  operator  the  effect 
of  an  occurrence  at  one  point  in  time  tends  to  bo 
spread  over  all  frequencies,  Furthermore,  time 
functions  with  corners  have  Fourier  transforms 
that  have  a  "ringing"  effect.  Thus  It  seems  de¬ 
sirable  to  have  a  data  window  of  a  high  continuity 


class  so  that  lh**  r^-Jifir’  sa;r.pl<-  function  •o»'‘ 
smoothly  t«>  /cm  at  it::  cuds. 

Such  cons  i  tut  at  i  r;otiv«<t.  eve  t  exa-ir** 
the  p s i b i  1  i  i  v  (>f  usi'v  (Jala  witVi-ws  based  upon 
higher  onhr  !*.-*. pi : -.e  basis  iu:.<  l  i  .  A  prelin- 
iiuit)  survey  is  tu.rou  i  «v  i  *•".  1‘iist,  we  i  ay  ciu»cs» 
the  dr  si  red  coot  i  uni  ty  class.  St  condl v,  1 1.  •  ctu  res¬ 
ponding  spevl ral  virduvs  have  th.  shape  ^  Hi — J 

so  that  the  magnitude  of  the  side  lobes  rapidly 
decreases  relative  to  the  rain  pi  ,d:  h'.  the  contin¬ 
uity  (or  order  k)  increases .  Tt  addition,  all 
spectral  windows  are  of  non -nev.a  t  i  ve  type  (Pfirxet* 
[13]).  fhirdly,  deb out  f  7  ]  lias  thrived  a  genetic  I 
n  1  gorithm  which  can  he  used  for  Accurately  evalu¬ 
ating  any  order  I>-spline  basis  function  with  nnv 
interval,  of  support  and  at.  any  nu  r.b  length. 

3.2  THE  ITNlVaMKlliAl.  SH  I*.:  ITNC.i  1  OSS 
Schoenberg  [  1A  ]  introduced  the  spline  funr- 
tions  in  19AA  with  the  folKwitiv'.  definition.  With 
k  a  positive  integer,  a  real  function  «A(l)  cu¬ 
ff  red  for  I  (.  b  is  called  a  spline  fuetiii’. 

of  order  k  or  degree  k-1  if  it  has  tK?  folic*- 
Ing  properties: 

(i)  A(t)  Is  of  class  c‘_!!(P), 

(if)  <7k(t)  fs  composed  of  polynomial  arcs  of 
degree  at  most  k-1  . 

The  polynomial  arcs  arc  Joined  nowhere  other 
than  possibly  at  Integers  n  If  k  is  even 
or  at  points  n  +  J  if  k  is  odd.  The 
spline  is  said  to  have  knots  at  these  points 
If  we  define 

,.k“i 


(HO 


yf1- 


i: 


y>  0 


y  <  0 


with  exception  of  Hie  special  case  0.  r  j  and 
let  6k  represent  the  k"'  older  central  diMT:- 
ence  operator  with  unity  step,  l.c.,  6~f(y)  * 

f(y)  and  gf(y)  -  6k-’f(y+?)  -  6‘"’f(y-'.) 
for  k  =  1,2,...,  then  the  function 

M*>  “  TiTTiT  Rkt* 

Is  easily  shown  to  he  a  simple  c -ample  of  a  spline 

function  of  order  k.  Mk(t)  Is  called  the  funda- 
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mental  B -spline  basis  function  of  ordt?r  k  .  Wo 
list  a  few  of  the  many  important  properties  of 
Mk(t)  and  refer  the  reader  to  Schoenberg  [r*,15], 
Curry  anti  Schoenberg  [6],  and  to  doBoor  [7)  for 


these  and  o' 


her  properties. 


k  k 

in  the  open  interval  (-  ) 


,  (t)  is  positive 
and. idont leal ly 


j  +  v  where  v 


zero  elsewhere;  it  has  knots  at 
is  an  integer  satisfying  0  <  v  <  k.  Mk(t)  is 
an  even  function,  and  its  integral  over  R  is 
unity.  Mk(t)  also  has  the  representation 


Mk(t)  * 


(k - 1) *  v=o 


k  k-V 

(-1) 


•h» 


r- 


For  example,  the  fourth  order  function  is  given 
by 

iti  <  i 


lU(t)  = 


i(2.  Ul)».  !)• 


,  t<  It  I  <2 


fotborwi se. 
We  shall  require  two  integral  propertied: 


and  -  MpM(t)  . 

3.3  HIE  PROPOSED  CLASS  OF  DATA  WINDOWS 
3.3.1  fill  roduct  ion. 

The  class  of  spline  data  windows  vc  propose 
consists  of  the  fundamental  spline  functions  each 
dilated  so  that  its  support  is  equal  to  the  length 
of  the  sample  function  and  normalized  so  that  the 
total  power  of  the  data  remains  invariant  when  the 
spline  data  windows  are  applied.  If  we  let  0k 
iepresc.it  the  normalization  factor,  then  we  can 
use  the  results  of  .the  previous  sections  to  ob¬ 
tain  the  following  table  of  cho  spline  data  win¬ 
dows  and  of  the  corresponding  frequency  windows, 
lag  windows,  and  spectral  windows. 

Since  the  total  power  >r  is  given  by 
nTf 

-•  *  x"(t)dt  *  T  Vnr[x|, 

‘T/r 

We  seek  constants  C*  depending  upon  the  order 
k  of  the  spline  data  window  such  that  also 


where  we  use  the  tilde  to  indicate  an  un -normal¬ 
ized  window.  The  oxprctai  ion  ot  this  equation  is 

<*  =--  c;  i'  ‘  EluquUKlx^dildt 

'-■■3 


=  C;;  v.u-(x|  ’  W;(t)<lt. 

‘-,/8 

Changing  the  variable  of  integration  we  obtain 

/  ^ 

c;  =  k/f  H?(t)dt  . 

>2 

An  application  of  Pars  oval's  theorem  yields 

•*  •  -AW 


Therefore,  we  can  calculate  the  appropriate  nor¬ 
malizing  factor  for  any  order  spline  data  window: 
Cl  “  1.0  and 


C. 


2  L-Alh hslllal!. 

.  i !  (2k-i-l) 


for  kS?. 


Norma  1  i7.  i  ng  for  invariance  of  total  power  is 
equivalent  to  the  requirement 

i  (Oi  :,JJ  Q(  f'df  -  i 

J  V' 


Data  Window 
W,(t) 

C,^M«(t)  -  <Vv(V) 

Frequency 

W  i  ndew 

1,(0 

C,  1'  /  .in(~ll7k)\k 

k'  \~n7k'  ) 

Equ  i  valent 

Lag  Window 
U(t) 

7-[k-  6v.MJ,(r)J 

v-(T) 

Spectra l 

W  i  nd  ow  ! 

%(i) 

(SK-mi- 

1 

r'i  /  >  hu  -n/k)  \  s‘ 

k;'  \  i/k  ) 

lAill  Jr.  1. 


[ cvff, ( c ) |  ot 

■  -V, 
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We  point  out  that  Wx  is  simply  a  trunca¬ 
tion  function,  f.j  is  the  lay;  window  suggested 
hy  Bartlett  ( l  ] ,  and  is  the  corresponding 

bare  Lott  spectral  window.  Wa,  I,a ,  and  are 
respectively  the  Parzen  data,  lap,,  and  spectral 
windows , 

For  the  problem  of  balancing  the  resolution 
and  stability  of  power  spectral  estimates,  ex¬ 
plicit  expressions  for  the  bias  and  variance  of 
smoothed  estimators  have  been  developed  in  the 
literaJturo.  Jenkins  and  Watts  1101  and  others 
give  these  expressions  for  the  indirect  computing 
method.  Welch  1 17, 18]  derives  the  expected  value 
and  variance  of  the  discrete,  directly  computed 
estimates  for  a  stochastic  process  which  is 
ilnussian,  ia  addition  to  our  previous  assumptions 
about  the  process.  In  each  case  some  property  of 
i be  window  being  used  is*  exp! ic i t ly  involved. 

Lh»  practitioner,  however,  needs  seme  device  to 
aid  in  the  select  ion  of  a  window  that  will  yield 
i  compromise  suitable  for  his  purposes.  The  most 
ce'.mon  measure  of  windows  is  c  tiled  bandwidth  (by 
malogy  with  filtering  problems)! 

Several  definitions  are  used  tor  bandwidth. 
Hackman  and  fukey  [1J  use  inst»ad  the  term 
equivalent  width,  which  they  define  to  be  the 
ratio  of  the  square  of  the  info  *rai  of  the  spec¬ 
tral  window  to  the  integral  of  its  square.  If 
we  define 

i  P  <i*(n<nt 

J  -30 

tlnn  the  equivalent  width  is  given  by 

b,  (<» 

tonkins  and  Watts  llOj  define  bandwidth  to  be 
the  width  of  the  rectangular  (speetral)  window 
for  which  the  variance  ot  the  indirectly  computed 
■'.p.*ttra1  estimate  is  the  same  .»•;  it  is  f-»r  the 
.'iven  window.  this  definition  of  bandwidth  %  i  «•  I  »1  - 
{  Ml) 

B,’<Q)  \  . 

As  has  been  noted,  tor  a  properly  normalized 
window  i<J(t)dt  1;  therefore  ^  .  We 


mention  Ilf  for  its  application  to  the  spline 
lag  windows  given  in  lablc  1.  However,  its  de¬ 
fining  piopeity  is  no  long*  i  title  for  tie'  direct 
method  of  computation  which  we  advocate.  For 
the  class  of  spline  windows  ve  have 


r.-i  (-n'a  - 

k'  ili-tk-i-l): 

thus  B ,  ( Q  •„  )  is  oas  1 1  y  detoraii  nod  by  ir.vorting 

tk  .  By  using  least-squares  techniques  tor 
k  =  3 , , . . ,  34 ,  we  oh  Coin  the  apjiroxiaiat  ion 
M'Q«)  .a  (l.f)S(k)u  0.32) /T  . 

Parzen  (13J  defines  bandwidth  .is  (ho  width  of  a 
rectangle  which  has  (.ho  same  area  and  sane  r.axl- 
iv.n.n  height  ns  the  S|'<  ot  r  >1  wirAw,  it.,. re¬ 
fute  ’ '  n  (i  t  d  r 

it.  io  i  — _ _ 

i-.inv  Q(f) 

and  for  the  'q.l  in,  wiu.h.ws  wo  havo 

O 

^  • 

Another  coin-only  ns«!  :;V*as>ir  ■  of  h  'advidt  h  is*  the 
distance  *vtwe«.n  the  ha  If -power  points  ot  the 
spectral  window,  for  an  oven  .:p..ciral  win¬ 

dow  contend  at  i  •  0,  N,(Q)  't  win  re  7  is 
the  f  i  e  ^  t  t  ri  quency  such  that  0<f)  1  Q  t0).  T.et 
j  n  l  •  i  »*«  ,  res-  :\t  rh,  *i:dtericiil  Hos.,ol  functions 
of  t  he  t i r«t  k i nd .  t  hen 


Approx  i-  it  rn  l>.  ;:-ch>  vi  i  h 

r-fe.e  ■  ,  t  '  taM.  i  ,  !  I  1  \  Also,  t 

qtien t  1  \  retiienc.  «  oe  :-ade  fo  ih.  .et  * i  i . >n 
hot  we -a  Mn  first  ;  ..j  -e».  .  t  v  » l  :  adovi . 

/lie  Zeroe-i  o|  »  «vCIT  it  t  t  -;hi«vp 

I 

fl  is  »  “i  itii!  f  I  4*  ’  •  h.-r . 

va»  i  >u  ■;  .uthois  h.i.-.  di  *«.n  thit  ike  ‘ins 
ot  sine- 'this!  - 1 .  ■  c  t  t  i 1  •  ■  f  i  *:ia  t  «  >  r  s  i;ur«  i  u*s  as  the 
bandwidth  n|  the  Wi  im-rrn  ..-s  ,  hut  that  the 

variance  Is  inv»-tsel.-  pioporlional  to  i  he  h.md- 
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.  dth.  It  can  also  be  argu  !  that  the  bandwidth 
t  ;ii.  wimtov  nisi  be  at  least  as  ssiili  as  the 
width  ‘-'I  the  nit  i  owe -r.t  i’itpj'M- at  detail  in  the 

»  t  I  Uvl. 

iuble  2  is  a  compilation  of  the  v.uieus  cen- 
*:t.iits  assooiatiJ  with  the  spline  windows  given 
x  it  i  a  :>*•■>  1  with  order  up  t  o  k-21  . 


Order 

Four  t ant 

band 

;idtl.*j  >  i 

k 

cu 

(QO 

•  <Q  .  • 

1 

1 .PT:0 

1.S00O 

1. 

? 

1 .  7  V/  0 

1.8S43 

1 .3313 

1.27S7 

2.3V.S 

2.301/ 

I . f f  dd 

: .  S7 !  6 

2.  f"  8  7 

2.  (>?/-/ 

1.91 V s 

I.R701 

8 

3.4PT3 

3.0203 

2  .  1  S2  1 

?  .  (1  <  •'  '» 

,  .  b 

3.9077 

3.3214 

2 .  i  .  3 1  ■ 

a  p  '}  ■,  y 

7 

4.3770 

3 .  S9 74 

2. s 474 

2.4 10s 

8 

4.8348 

3. 8  Sit 

2.7170 

?. .  Srf>3 

0 

S.2V83 

4.0930 

2. O071 

7  '*Vt 

M 

s  .  7  1  07 

A.31HW 

3 . 0669 

2.80 10 

11 

(>.131? 

/,  .  S3t»3 

3.2188 

1 . 0  1  ‘  1 

17 

6.S42  7 

4.7423 

3.31.3** 

3  .l/oi 

n 

6,9*  ’:>$ 

4,9394 

3.4030 

3.3011 

14 

7.3'.  13 

S . 1292 

3.4387 

3.4264 

IS 

7.7298 

S. 3121 

3.7087 

3.  S/,73 

16 

8.1118 

S.489  0 

3.8904 

3 . 604 1 

17 

8.4880 

8.6602 

4.0113 

3.7774 

18 

8.8S87 

8.8268 

4.1286 

3.8873 

19 

9.2243 

S.Q882 

4.2427 

3 . 9>'4  3 

20 

9 .  SS SI 

6 .  14  Sf, 

4.3S38 

ft,  00  8  i 

21 

9, 9'.  14 

(-.2900 

1* .  462 1 

— 

■t .  7  O00 

TABIE  ?  . 

3.4  CONiTTINC  WITH  SDL1NE  DATA  WINDOWS 
Except  foi  the  choice  of  data  windows,  the 
method  we  propose  for  computing  power  spectra  is 
that  of  Welch  (18).  Departing  from  our  conven¬ 
tion,  we  assume  that  we  now  have  a  4 ample  func¬ 
tion  x(t)  wr  the  interval  0<t^S.  Using 
for  example  a  bandwidth  criterion,  choose  the 
order  k  of  the  spline  window  to  he  used  and 
then  the  vfnd'«’  length  T.  Subdivide  y(t)  into 


N  segment's  each  of  length  T  and  with  a  dis- 
pI.v«iaou  D  between  the  initial  points  of  n«!i- 
«e<  til  j  ve  seg::«iils.  Denot  c  t  lit'  S\  y  Milts  by  >; ,  ( t  )  . 


Thus  (N  -  I)!;  t  V  -  S,  and 

x.(t)  * I  (j  -  !/!>8  lor  j  I . K;  f;  . 

Welch  nsi  s  the  clti'lce  of  I)  unlv  as  a  parameter 
to  increase  the  '.rd-ilitv  of  the  resulting  esti¬ 
mated  power  spec  Li  urn,  In  the  c.K»  ol  split*..  win- 
do..*s  ;m  additional  consideration  should  he  made. 
11  I)  -  T ,  some  oi  the  data  is  weighted  substan¬ 
tially  heavier  than  the  remainder .  As  c an  be 
Seen  from  Figure  1,  the  inequity  is  more  pro- 
no.. need  as  the  order  of  the  spline  window  in¬ 
creases  in  the  sense  that  the  spline  data  win¬ 
dows  become  progressively  more  peaked.  Since 
oin-  normally  assess  that  all  portions  of  his 


data  aie  equally  representative  of  tin  stochastic 


process,  intuitively  the  order  of  the  spline 
window  should  also  influence  I>,  In  practice  a 
good  rule  of  thivb  for  relatively  short  s.v.ple 
functions  seems  to  b^  D  -  -  T  for  k  "  A.  I)  ^  T 


1  c' 
foi  k  -  10,  and  1)  ^  7 T  for  k-  ,':0  with  t he 

ovei  lap  for  other  order  window:'  t.i/.id  on  those  . 

c 

Next  compute  the  smoothed  spectral  estimates 

r.CO  ’  x,(t)w.  (t)o'2nU,dt !“,  j-1 . N. 

1 -'/*  ■ 


Finally,  the  estimated  power  spectrum  is  given  by 

r.(f)  -  i  >■  P)(n  • 

w  i  - 1 


Regardless  of  the  order  of  spline  window 
used,  P,(f)  is  obviously  a  consistent  estima¬ 
tor  of  P(f)  (at  least)  whenever  the  stochastic 
process  is  Gaussian,  'lo  see  this,  the  express¬ 
ions  for  E[P#( 1)1  and  Var[re(f)]  as  given  by 
Welch  117,18]  need  only  be  compared  with  the 
convergence  properties  of  Lhe  spline  windows  as 
discussed  in  this  paper.  That  i ,  as  S  s*  ve 
have  the  following  type  convergence  for  a  sta¬ 
tionary  Gaussian  process:  if  T— a  such  that 

£-0,  then  E(  <P<  f>  -  F,<r))a)-0  . 

In  implementing  this  procedure  on  a  digital 
computer,  one  could  derive  a  closed  form  for  the 
R-'-jlfne  basis  functions  like  that  given  pre- 
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cov.irinnci*  fund  ion. 


viously  for  However,  for  .1  general  p^o- 

grfln  am!  one  that  won  hi  undoubtedly  bo  numeri¬ 
cally  preferable,  the  algorithm  of  deftoor  [7]  is 
roeiw»mended  for  evaluating  the  spline  functions. 

Jt  couts.',  by  symmetry  only  half  the  window  need 
be  directly  evaluated.  Brigham  (''♦I  thoroughly 
discusses  the  use  of  the  K7-r  for  output  ing  the 
discrete  Fourier  transform  which  approximates 
the  Fourier  integral  we  used  in  writing  Pj(f). 

4 .  EX  PE  R  l  MENTA  T.  RES  U  LTS 

The  effectiveness  of  the  spline  data  windows 
is  readily  apparent  from  experiments  with  machine 
generated  autoregressive  processes,  power  spec¬ 
tral  estimates  computed  using  various  order  spline 
data  windows  are  compared  with  those  which  used 
the  Bingham -God  froy-Tiikey  and  Parson  data  windows 
and  with  estimates  computed  by  the  indirect  method 
of  dlnokrun  and  Tukov  incorporating  the  Hanning 
iag  windev  [31.  The  estimates  were  not u..i  1  i/.ed  bv 
dividing  out  the  viri.ece  of  !he  to«*h  .sMc  pro¬ 
cess;  as  is  customary  when  one-sidid  power  spectra 
a**e  used,  the  magnitude  has  been  doubled. 

For  these  experiments  we  used  sample  func¬ 
tions  of  -nitoregress Ive  processes  which  were 
machine  generated  according  to 
3 

xt  Oj  *.  - .  I  «t  ,  t  -s  1 , . . . ,  Ltj24 

to 

where  m  is  the  order  of  the  process,  f«>t }  are 
the  regression  coefficients,  and  is  a  white 

noise  pioCi-ss  with  E{-\  |  U  and  Vnr[s,|  -  1. 
rhe  parti«;ilnr  process**;  which  were  used  in  the 
i  I  lust  r  ir  ioiv:  had  the  following  coefficients: 

(i)  i-’icsf  order,  -O.nO 

(ii)  Second  order,  »*,  1.00,  -0.50 

(iii)  Fourth  order,  0.73,  » ,  0.30, 

r,  '  -0.33.  ...  '  -0.30  . 

!he  following  scheme  identifies  the  windows  in- 
v«j!  /••<!  In  » tie  c -.•input  nf.  ion  of  the  virf.mi  power 
spectral  *•  *  t  i  mates  . 

(i)  dinnin.c  .\x:  estimate  computed  by  the  in¬ 
direct  s.i  ihwd  with  the  Hanning  lag  window 
[}J  and  using  x:<  lags  of  the  sanple  ;mto- 


(ill  h-G--  King!,  a.  1  t.od  f  re/-'!  wke  d.,i  a 

window  was  use*!  vi  tfi  segments  of  the 
:»a.:.*)Je  function  of  length  T  v;< 
samp  1 es . 

(lli)  Parboil  ss:  Paraon  data  window  was  used 
Willi  so psion is  (it  tho  sa.:.plo  function  of 
lcu:;t!t  T  ■-  xx  simples. 

(iv)  tip  1  ino  xx-yy-K<s:  Spline  data  window  of 
order  k  =  xx  was  used  with  so)>r.r  :it  s  of 
••ii.-.ple  fund. ion  each  with  length 
T  J  yy  samp  It's  and  d  i  sp  I  acotaont  D  ■■  tt 
sarnies  between  initial  points, 
iho  graphs  i t ru 1 1  r.e  these  experiments  arc 
found  in  Appendix  1, 
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